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ABSTRACT 


The equations for the kinetic modeling of a supersonically flowing 
electrically excited laser system are presented. The work focuses on 
the use of diatomic gases, in particular carbon monoxide mixtures. The 
equations presented include the vibrational rate equation - which describes 
the vibrational population distribution, the electron, ion and electronic 
level rate equations, the gasdynamic equations for an ionized gas in the 
presence of an applied electric field, and the free- electron Boltzmann 
equation including flow and gradient coupling terms. The model developed 
accounts for vibration-vibration collisions, vibration-translation collis- 
ions, electron-molecule inelastic excitation and superelastic de-excita- 
tion collisions, charged particle collisions, ionization and three body 
recombination collisions, elastic collisions, and radiativq' decay, all 
of which take place in such a system. A simplified form of the free- 
electron Boltzmann equation is developed and discussed with emphasis placed 
on its coupling with the supersonic flow. A brief description of a 
possible solution procedure for the set of coupled equations is then 
discussed. 


INTRODUCTION 


The concept of a fast flow laser system is not new . It has been 
demonstrated for electric discharge C0„ laser mixtures, in both cw and 
pulsed operation, that a considerable increase in power can be achieved 
with the use of fast flow 2-11 . The addition of fast flow gives several 
attractive features! (a) the power output scales with the mass flow 
rate - this can represent up to a factor of 10 5 increase in power over 
that obtainable from conventional diffusion controlled laser systems, 
(b) the translational temperature of the gas is maintained at a low 
level by the convection which prevents the lower laser level for CC >2 
systems from becoming populated to such ani extent that the inversion is 
destroyed, and is very important for diatomic laser systems, which we 
shall mention below. 


^Research Scientist. 

The author would like to acknowledge R.L. McKenzie and Dr. J.R. Viegas 
of Ames for valuable discussions and information. This work was carried 
out under National Academy of Sciences — National Research Council 
Research Associateship . 



To understand the importance of low translational temperatures for 
diatomic laser systems the mechanisms responsible for diatomic laser 
action will briefly be mentioned. When considering diatomic molecules, 
the anharmonic structure, i.e. the small decrease in spacing between 
vibrational energy states with increasing quantum number, introduces a 
new and crucial factor that leads to the processes which can ultimately 
produce laser action. Treanor et.al. have shown that during the 
collision of two such anharmonic diatomic molecules the subsequent exchange 
of vibrational energy (V-V exchange) can lead to significant departures 
from a Boltzman distribution of vibrational energy. Specifically, in 
Ref. 12 it has been shown that under the conditions of a low translational 
temperature together with a high degree of vibrational excitation, 
population inversions between the combined vibration-rotation levels of 
adjacent vibrational states can occur. Such inversions then provide 
the basic mechanism by which diatomic lasers operate. In addition, to. 
maximize tt*e radiative gain on a particular vibration-rotation transition, 
once the inversion has been achieved , the number of molecules at the 
level in question must be maximized. This occurs when the molecules are 
preferentially populating the lower rotatiohal levels, which occurs for 
low translational temperatures (assuming rotation is in equilibrium 
with translation) . 

The conditions of high vibrational excitation and low translational 
temperature occur naturally in the supersonic expansion of a gas from 
a hot reservoir. The occurance of these conditions in such gasdynamic 
expansions led McKenzie to develop the CO gasdynamic laser. Defining 
efficiency, for the gasdynamic laser, as the -ratio of laser power output 
to input enthalpy flux, where the latter is the power used to produce the 
heated reservoir and associated supersonic flow* McKenzie obtained 
efficiencies on the order of a few tenths of a percent for the CO system. 

On the other hand, electric discharge lasers without flow have very 
high efficiencies but low power. For a CO system Bhaumik. et.al. has 
recently obtained 46% efficiency for conversion of electric power to 
laser power. This is the highest efficiency demonstrated for any laser 
system, polyatomic or diatomic. For such a system the laser gas is 
essentially in a static condition, cooled by liquid nitrogen, and excited 
by a glow discharge. The reason such high efficiencies are possible is 
that the input electrical power goes predominately into exciting the 
vibrational levels of the gas by electron-molecule collisions. This 
leads to the condition of high vibrational excitation, while the liquid 
nitrogen cooling provides the low translational temperature. However, 
since this is a static discharge laser, and consequently diffusion 
controlled - i.e. removal of heated gas is by diffusion to the walls, the 
advantages of increased power levels with fast flow rates are not utilize^. 

Consideration of these facts indicates that the use of electric 
discharge excitation combined with fast flow rates could potentially lead 
to highly efficient high power lasers. This is based on the separate 
advantages of each, and assuming that when combined they do not adversely 
effect each other. 
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Recently, laser power, using two different schemes but employing the 
features of electric discharge excitation and fast flow in CO mixtures has 
been demonstrated. First, Rich et.al.-^-^ (and later Kan et.al. °) excited 
a CO mixture with a glow discharge in the plenum chamber of a supersonic 
nozzle and then allowed the excited mixture to expand to supersonic 
speeds. It is clear that such a scheme produces the necessary conditions 
mentioned before for diatomic laser operation. However, this scheme 
produced a relatively inefficient laser with efficiencies of the order of 
a percent. This is due to the fact that much of the vibrational excitation 
produced in the plenum is lost to heating the gas in the subsonic portion 
of the expansion. The second scheme, developed by Kan et.al.l?, is 
essentially a direct extension of the liquid nitrogen cooled static 
electric discharge mentioned previously, whereby a fast flow rate of 
CO is now employed. This scheme used subsonic flow rates and relied 
on liquid nitrogen cooling and thus did not utilize a supersonic expansion 
to give the low temperature, and the concurrent further increase in 
power that would result from such an increased speed. However, this 
scheme did achieve efficiencies up to 26%. 

These two schemes, representing initial approaches to combining fast 
flow and electric discharge excitation in diatomic gases have certain 
disadvantages, as mentioned above. With the prospect of more fully 
utilizing the advantages of each feature the intent of this work will 
be to discuss the situation of electric discharge excitation occuring in 
a supersonic stream. 

This leads to important questions as to the mechanisms taking place 
in such combined systems, and the result of using supersonic flow, 
with the discharge occuring in the flow. For the case of electric 
discharges in a supersonic stream, questions concerning the stability 
must also be answered. Here, we assume a priori that such a discharge 
is feasible. 

With these ideas in mind it is the purpose of this note to develop 
a kinetic model for the case of a supersonically flowing CO mixture 
with a stable electric discharge established in the flow. This note is 
in the form of a summary of the work done to date which consists of a 
development of the equations relevant to the problem and a brief 
description of a possible solution procedure. 


BASIC THEORETICAL MODEL 


Physical Model 

The geometry employed is that of a supersonic nozzle in which the flow 
is to be established (see Fig. 1). In some region of the nozzle, electrodes 
are situated such that a discharge, alligned with the flow can be produced. 
The laser power is then extracted transverse to the flow and discharge 
direction. In the sketch shown, the discharge is located in a region of 
uniform flow properties, although it may be located in any portion of the 
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nozzle; even the reservoir, as in Rich's 15 experiment. 

For efficient excitation of the vibrational levels of a gas by electron 
collisions the average electron energies are those conveniently obtained 
in the positive column of a normal glow discharge (around 1 ev) . We shall, 
therefore, consider only glow discharges, assumed to be fully established 
in the flow, and then attempt to model the ensuing positive column region. 
The questions of instability of the discharge are not considered. Such 
questions are, of course, important but can be investigated separately. 

(See Ecker et.al. 18 and Nighan et.al. 19 ). By limiting the modeling to 
the positive column the space charge effects at the electrodes can be 
neglected. They are crucial to the establishment of a stable glow 
discharge, but unimportant in the physics of the positive column. . The 
problem will also be assumed to be convectively controlled (see Section 3); 
thus the effects of diffusioi^to the walls will be considered negligible 
and the standard Schottky 20 ’ 1 type analyses are not relevant. 

The gas is further assumed to enter the discharge region in a pre- 
ionized state. Preionization has an important effect on efficient laser 
operation. Ordinarily, without preionization, the glow discharge is 
required to be self-sustaining. Then, ionization, produced by electron- 
molecule collisions at a rate sufficient to balance the loss mechanisms 
present, takes place within the discharge itself to keep the gas con- 
ductive. However, schemes that separate the above function of ionization 
from the concurrent vibrational excitation (crucial to. producing the 
inversion needed for lasing) have been shown22»23 to considerably increase 
laser performance. This is because the electrons in the discharge that 
are produced by an external means, can be of lower energy than that 
required to sustain the discharge. In fact, by proper tailoring of the 
discharge voltage, the electrons in the discharge can then be used for 
efficient vibrational excitation exclusively. 

The emphasis of this work is thus to model the interaction of a 
(preionized) supersonic flow with the plasma of a positive column to 
study its effects on the resulting vibrational population distributions. 

Energy Transfer Mechanisms 

In a discharge of the type considered, several collisional and radia^- 
tive energy exchange mechanisms take place. These transfer mechanisms 
are conveniently modeled by a division into six energy modes, consisting 
of a free-electron mode, an anharmonic oscillator mode of neutral 
particles, a bound-electron mode of neutral particles, a kinetic mode of 
heavy particles, and a UV and IR-photon mode (see Fig. 2). Details of 
the resulting model are similar in some parts, to the models described 
in Ref. 24, 25, 26. 

Referring to Fig. 2 we assume that the free electrons, accelerated by 
the applied electric field, participate in the (numbered) energy transfer 
mechanisms of: (1) inelastic excitation and superelastic de-excitation 

collisions with the vibrational and electronic levels, (2) elastic 
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collisions with heavy particles, (3) ionization and three-body recombination 
collisions, and (4) electron-electron collisions. This energy distribution 
is subsequently redistributed by: (5) vibration-translation collisions, 

(6) vibration-vibration collisions, (7) radiative decay of electronic 
levels (8) radiative decay of vibrational levels and finally, (9) possible 
induced radiative transitions between vibrational levels. Several 
assumptions are implicit in this model: (a) the gas mixture consists of 

electrons, singly charged ions, and nuetrals - of which several monatomic 
and diatomic species are allowed, (b) translation and rotation are in 
equilibrium at the temperature T, (c) the excitation of ion electronic 
states is negligible, (d) the ions are in kinetic equilibrium with the 
neutral particles, (e) dissociation is ignored, and (f) photoionization 
and two-body recombination are negligible. 

One further simplification will be to neglect at the outset the effect 
of stimulated emission (process 9) since we are only interested here in 
the calculation of radiative gain prior to the introduction of any laser 
applications. 

The free electrons in this model are not expected to be in a Boltzmann 
distribution. In fact, for conditions typical of CO and CO^ lasers, it 
has been shown theoretically^ and experimentally^’ , the distribution 
is profoundly non-Maxwellian and is strongly affected by the electron- 
molecule energy processes discussed above. Since the electron-molecule 
collision rates (needed to determine the vibrational population distribution) , 
the electron temperature, the heat conduction, and the current density, are 
all averages over the distribution of electron energy it is clear that a 
knowledge of this distribution is needed for an accurate representation 
of the model. 

i' ; 1 

In summary, the energy transfer mechanisms discussed above, together 
with an appropriate development of the f ree-electron-Boltzmann equation 
are postulated to sufficiently model the combined flowing gas-electric 
discharge situation under consideration. In Section 3 to follow, the 
necessary rate equations for the model are given. Section 4 presents the 
standard gasdynamic equations • for an ionized gas in the presence of an 
applied electric field. Finally, in Section 5 the free-electron Boltzmann 
equation is developed and certain simplifications are then made that are 
appropriate to the present model. 

RATE EQUATIONS 

The equation necessary to describe the vibrational population dis- 
tribution of the laser gas mixture is the species continuity^, 31 e q ua ti on 
for the vibrational levels. For a general unsteady flow problem this 

equation is '■ 

. k 

' -gj* + V • (n* v s ) = n£ k = 0,1,2 ... N s , (1) 

1c 

where n g is the number density of species s in level k, v g is the 

species velocity and, h k is the net rate of production of species s 
in level k to be determined from a consideration of the energy transfer 
processes described previously. 
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• k 

Before proceeding with the expression for n the conditions that 
insure a convectively controlled problem will be presented. We find, 
after an examination of the left-hand side of Eq. (1), that two conditions 
must be satisfied. These are: 


and 


u 0 » 


_L 

d 2 


D 


amb 5 


( 2 ) 


u Q » v f , 

where u , V , L, d, and D ^ are: characteristic flow velocity, electrpn 
diffusion velocity in the a supersonic flow direction, discharge axial * 
length, discharge transverse length, and ambipolar diffusion coefficient . 

Each condition can easily be checked once a solution is known. In 
fact, values for V^, can be calculated exactly from the formulation to 
be discussed in Section 5. _ 


For the present problem we assume that the above conditions are 
satisfied and thus Eq. (1) can be written for each vibrational level k, 
for steady one-dimensional flow, as follows: 


P u 


dy \ p ) 


k = 0,1,2 


No 


(4) 


where y is the coordinate in the one-dimensional flow direction, p 
is the average density, and u is the mass average velocity. 


To determine the expression for the production of excited vibrational 
levels, n , we take into account certain of the energy transfer processes 
described S in Fig. 2. These consist of the following collisional and 
radiative processes: 

(i) Inelastic and superelastic collisions with electrons - process 

(ii) Vibration-translation collisional exchange with the various , 
constituents - process 5. 

"(iii) Vibration-vibration collisional exchange within the laser gas 
and other diatomic constituents - process 6. 

(iv) Excitation due to the spontaneous emission from the excited 
electronic level - process 7. We assume here that only one 
electronic level (a 3 II) is excited sufficiently to participate 
in the energy transfer processes. 
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The development of these conditions follows a discussion given by 
Hassan and Bordeaux . Here, however, the diffusion to the walls takes 
place by electron-ion ambipolar diffusion-^. 


(v) Spontaneous emission - process 8. 

(vi) Possible ionization and 3-body recombination effects. 
In addition, the modeling includes the following features: 


(vii) The vibrational quantum states of the diatomic gas, are 
treated as those of a Morse anharmonic oscillator"' : 


e k = K0 s k[l-A s (k+l)l , 


( 5 ) 


where 


9 = characteristic vibrational temperature, 

s 

A = small anharmonic oscillator parameter, 
s 

k = Boltzmann constant, 

Ic 

= energy of the k vibrational state of species s. 

(viii) With the exception of vibrational transitions induced by 
electron collisions, only single quantum transitions are 
allowed. 

(ix) The free electrons in the discharge are not assumed to be in 
a Maxwell-Boltzmann distribution and consequently several 
parameters will be needed to characterize their properties. 

(x) For the excited electronic state its ground vibrational level 
is the only level assumed to be populated. 


Combining these various features of the model with Eqs. (4) and 
defining the mole fraction of species s in level k as 

,k _ _£i_ 


- y k » 

Zn s 

s ,k k 

we obtain the following vibrational rate equations for x : 


( 6 ) 


dx| 

dy 


= £ Z sp Q sp + x e | £ x s 


k- 1 


e ^ek 


k-l 


8=0 
N c 


8 = k + 1 


x® Rs k 


-x s k £ R ks -X k £ 


R kS - X S K Rs A + xj Xe R; 


>k \ 


C s /\ e 


8 = 0 


- x? A 


k,k-l 


+ X 


8 = k+ 1 
k+1 a k + 1 ,k 


+ xf A 


Ek 


k = 0,1 ,2 . . . N s 


(7) 
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In this equation the first term on the right-hand side accounts for vibration- 
vibration (V-V) and vibration-translatio^ (V-T) collisions and is given, 
between species s and p, by McKenzie as follows: 


Z 

sp 


average 

= P 
m 


collision rate 
a 2 I 8ttkT 

dsp v "ST.; 


( 8 ) 


where d 


2 

sp 


average collision cross section for encounters between 
species s and p. 


u = reduced mass , 
sp 


m = molecular mass. 



Im , 
s s 


molar production per collision 


= [Qsplv-V + [Qsplv-T ’ 


( 9 ) 


( 10 ) 


where 


and 



(ID 


( 12 ) 


In the above notation P 3,b is the V-T probability per collision that 
species s will make till transition from vibrational quantum state a to 
b as a result of an encounter with species p. Similarly P gp ( a ’ b ■) is 

the V-V probability per collision that an encounter between species s 

and p will cause an a to b transition in species s and a simultaneous 
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•V 


in to n transition in species p. 

It is beyond the scope of this paper to discuss in detail the appro- 
priate expressions for the above transition probabilities and consequently 
we will briefly describe those to'be used and cite references for the 
interested reader. For the V-T probabilities, which are reasonably well^ 
established, the functional form given by Schwarz, Slawsky, and Herzfeld 
modified by Keck and Carrier 3 ^ and normalized to fit the data of Millikan 
and White 3 ' are to be used. The expressions obtained in this manner can be 
found in McKenzie and also Rich et.al. The V-V probabilities to be used 
include contributions from a short range interaction as done by Rapp- 
Englander-Golden 3 ® and a long range interaction important at low temperature 
as treated by Sharma and Brau 3 ^. A detailed discussion of the expressions 
obtained from such considerations is given by Jeffers and Kelley . and also 
by Rockwood 


For the remaining terms in Eq. (7) the following definitions apply: 


x = mole fraction of electrons, \ 

e 

x^ = mole fraction of ions of species s, 

s ... 

E 

x = mole fraction of neutrals of species s in the ground 
s vibrational level of the excited electronic state, where the 
excited electronic state is denoted by E, 

ah —1 

A = Einstein A coefficient (sec ) for spontaneous emission for 

s species s, obtained from Penher- - for a Morse oscillator 



and is given by: ' 

» k.k-1 _ 

'ei-ef-V'. 

|R k,k 1 1 

2 ' ' ' 

A 10 




< e' / 

U;° 1 


where 

R k,kl _ / 1 2 \ 

r oh i 

1/2 (-l) 2k -> 

ho- 

2 k- 1 ) (i-2k-3 y 


R s 10 \ A < / 

L(H (H_ 

C'i 

*1 . 

H 5 


OH J 


(13) 


J (14) 

-J ' ” * / — 

•pt- — J. r 

A = Rate coefficient (sec ) for spontaneous emission from the 
s ground vibrational level of the excited electronic state 
to vibrational level 
species s^ 3 . 


64ir 4 


gk 3h 4 c 3 


(e kE ) 3 I Res I 2 q° k 


k of the ground electronic state of 

(15) 


electron-molecule vibrational excitation rate (sec ^) 

= £ f ef°(e,y)S^\e)de , 

m As 


(16) 
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R = inverse electron-molecule de-excitation rate 
s „ 


L 


eS k6 (e)f°(e-ef,y)de, 


R = collisional ionization rate coefficient 
S no 


= p / ay 


ef°(e,y)S kx (e)de, 


where 


R = rate coefficient for inverse three-body recombination process 

s ■■ - 
<■0 kX. 

J 2 f°(e-e^,y) 


= Frank-Condon factor for transition from ground vibrational 
level of excited electronic state to vibrational level k 
of ground electronic state, 

] Re g | 2 = electronic transition moment, 

kE 

: = energy difference from vibrational level k of ground 

“ electronic state to ground vibrational state of excited 
electronic state, 

e = free-electron energy (ev) , 

f°(€,y) = energy and spatially dependent isotropic electron 

distribution function (to be discussed in Section 5) , 

lcJi 

S (e) = energy-dependent collisional cross section for vibrational 
S excitation of species s, 

kX 

S (e) = energy-dependent collisional ionization cross section 
of species s , 

= e 8 -e k 

c s c s c s > 

e k \ = _ e k 
C S C S C S 5 


where e is the ionization energy, 
s 

In the above expressions detailed balance has been used to relate the 
inverse processes^. The spectroscopic constants needed can be found in 
Herzberg^ or in the NBS^^ tables. 

We note that if the electron distribution function is Maxwellian at a 
temperature T g then the inverse rates arfe simply related as follows : 
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( 20 ) 


and 



( 21 ) 


where the barred rate coefficients are calculated for a Maxwellian dis- 
tribution. 


Equation (7) is written for diatomic gases of species s. If a 
monatomic gas is also employed in the mixture then an appropriate rate 
equation, taking into account electronic excitation, ionization, and 
possibly spontaneous emission is required. This is easily obtained from 
Eq. (7) by setting the first and last terms on the right-hand-side to 
zero and interpreting the indicies of the remaining terms as associated 
with electronic levels of the monatomic gas.! It is clear, however, that 
with the magnitude of the electron energies that are dealt with in this 
model only one or two electronic;: levels will- be excited, with the pre- 
ponderance of the atoms being in. their ground electronic state. 

In addition to the vibrational rate equation, rate equations for the 
excited electronic level, the ions, and the electrons are also required. For 
the present model these equations are: 


electronic level rate 


dxf 


'-‘'’■s E nEo , /•y k\noE V E pEk . „2 I pkE E v aEE 

u = -XjXj Rs + X e (2. X s ) K s -X e X s K 5 +X e X s K s x e x s ^ , 

ay k i K 


( 22 ) 


ion rate 
.1 


EpE\ _ v 2 v 1 r \E = ^1 > 


dx s _ y k nkk v 2 Y 1 T + X R tA - — X 1 R 

U — “ X e L X s Kj- X e Xg Zj ix s « Xg Xg Ixg Ag Agi\g 

ay k k 

electron rate 

dx e v , • K 

u d7"-? (w * h 

where 


(23) 


(24) 


r Eo = rate coefficient for de-excitation by electron collision to 
s 1 'ground electronic state-assumed independent of vibrational, 
states. j )■ 


GASDYNAMIC EQUATIONS 

To determine the flow and thermodynamic properties the overall mass, 
momentum, energy, and the equation of state are needed. For the reader 
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who may be interested in problems other than the steady quasi-one-dimensional 
model employed here, the mass, momentum, and energy equations will first 
be written in general form and then specialized to the present model. 

Assuming an inviscid flow in the presence of an electric field with no 
imposed magnetic field these equations are^9» : 


overall mass 

3p 

at 


+ V • (p u ) = 0 , 


(25) 


overall momentum 
Du 


Dt 


+ Vp = p e E , 


(26) 


overall energy 


D / 1 2 v, 

p Dtl2 U +h 


_|P = -V • q + E- J + Q,n 


(27) 


where 

D _ 3 , rr 

Dt at u v ’ 

p = static pressure of mixture, 

P e = free charge density = (En^-n^e, e electron charge magnitude, 
h = total enthalpy of mixture, 
q = heat conduction, 

J = total current density = p e u + e(2 n* V* - n e V e ) , 

= ion diffusion velocity, 
s 

V = electron diffusion velocity, 
e 

Q = heat added per unit volume due to inelastic effects, 

E = electric field strength. 

We' now specialize these equations to steady quasi-one-dimensional flow, 
add the equation of state, and make the additional assumptions of (i) 
neutral plasma, (ii) mixture of perfect gases, and, (iii) no heavy 
particle heat conduction, to obtain the following: 


overall mass 

d_ 

dx 


(pAu) = 0 , 


(28) 
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overall momentum 


du x dp A 

ST + -fit = 0 • 


(29) 


overall energy 


pu dx 4 u2 + h) = + EJ + ^ 1N ’ 


(30) 


state 


-M 1 


P = ± {x e f£T e + KT [x e + 2 (xf + 2xf) + 2 (2xf)] 

s=d k s = m k 


(31) 


h = - + 4- l4x e »cTe + |x e KT+ 2 [x e *T + ^kT (2x k ) + 2x k e k + xf ef ] 

/ ^ s=d ^ 


m 


+ 2 [ 4 ifT (2 Xs ) + 2 xf ef ] \ , 

s=m * k k 


(32) 


where the nomenclature used Is: 
d = diatomic constituents, 
m = monatomic constituents. 


A ■ A(x) cross sectional area, 
k 


e = electronic energy of the kth level of monatomic species s, 

s =■ m 1 


and 


W N S -1 N s 

E E i« k -*^R 

k e=k+l 


f S + xf A‘ k ]ef s 


+ [(2xf ) R s oe - xf Rf° ] e° E + 2 [- xf R s kk + x e x'Rf k ] ef 

k k 


(33) 


+ [-xfR s Ex + x e x‘Rf E ]efM + 2 2(Z sp Qf p ef e ) 

J k p 

k*k-l 


To complete the formulation expressions for the current density J, 
electron heat conduction q , and electron temperature T are needed. 

The current density for a neutral plasma is given by the conduction 
current density as 

J = e(2n!vl-n e V e ). (34) 
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But for the case where the electron energy is much greater than. that of the 
heavy particles, such as is found in glow discharges V /V <<|, therefore 

the current is carried by the electrons only and is given by 

J » - en e V e . (35) 


The electron diffusion velocity V is calculated from an integral of the 
electron distribution function as 6 follows^: 


V e 


1 

3 



J ef'(e,y)de. 
o 


(36) 


The same is true for 


T^ and and they are given by : 


T e 

and 

Qe 


o° 

J I e 3/2 f°(e,y)de , 
o 

^ ~~ — L = x e f e 2 f'(e,y)de. 
3 m y/m e J Q 


(37) 


(38) 


In these expressions, as well as the rate coefficients discussed in the 
previous section (Eqs. 16-19) , the functions f (e»y) and f (e,y) appear. 
These functions represent the first two terms of a series expansion 
of the free-electron velocity distribution function. (The procedure by 
which they are obtained, and the equations governing them will be given 
in the next section.) They are obtained from a solution of the free- 
electron Boltzmann equation which is itself coupled to the gasdynamic 
equations . 


For the present discharge model we will be concerned with weakly 
ionized plasmas, that is t 

x e /2 Xs << 1 • ' 

s,k 

For this condition we find, after a comparison of the various cross 
sections that appear in the rate equations, that all three-body recom- 
bination terms can be neglected. This, coupled with the assumption of 
a neutral plasma eliminates the necessity of determining the mole fraction 
of ions so that the ion rate equation can be dropped from the system of 
equations. 

We thus have a system of 8 + s(l+N ) coupled equations consisting of 
Eqs. (7), (22), (24), and (29) through S (32) , and the equations for f° (e,y) 
and f 1 (e,y) (see Eqs. (42) and (43) of the next section) together wit^ 
auxiliary equations, for the system of unknowns: p, u, T, p, h, x e , Xg , 

Xg, f° (e,y) and f (e,y). 


Finally, when this system of equations is solved and with the solution 
for the vibrational population distribution known, the small signal gain 
can be calculated. A detailed derivation of the necessary gain equation 
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for the case of Doppler and pressure broadened lines for a diatomic gas 
is given by McKenzie* and will not be repeated here. 

ELECTRON KINETICS 


Boltzmann Equation 


As previously mentioned the free electrons in the discharge are 
generally in a non-Maxwellian distribution. The fundamental equation 
governing this distribution is the Boltzmann equation. Analysis of 
various features of this equation have been done in the past » » ’ 

however, it is convenient here to summarize those features appropriate 
to the present model. 


In a frame of reference moving with thq mass average velocity of the 
fluid, the electron velocity distribution function f(c,x,t) satisfies 
the equation* 


^(n e O + c" • V x (n e f) - 



n e V c f-c 'V x u-n e V c f 


df I 

'COLL ’ (39) 


where V and V are the gradients with respect to the spatial 
coordinates and eSectron velocity, respectively^ ^The right hand side 
of Eq. (39) represents the rate pf change of f(c,x,t) in a series as 
follows: 

f(c ;x- 1) = f°(c,x ,t) + f • f (C,x,t) + ... , ( 40 > 

where 

|T'(c,M)|«f 0 (c,x,t). (41) 

The function f° represents the isotropic contribution, and •?’ gives - 
the small anisotropic contribution to the velocity distribution function. 


Substitution of Eq. (40) into Eq. (39), and spgcessive integration 
of the velocity moments of^the resulting equation , yields a sequence of 
coupled equations for f°, f'. . . • • Transforming this set of equations 
to the independent variable e, the electron energy e= 1/2 m^c ; using 

the condition of Eq^ (41) and assuming that |u|fc, we obtain two coupled 
equations for f°(e,x,t) and f*(£,£,t) as follows: 


rem e \' /j D . ,o, , e „ , ~t., n e eE 3_ f f ■ ) 
-Y) Dt (nef ) 3 Vx(n e f ) - 3 3e (et } 


('■?) 

-(?) 


Du 

Dt 


^(e f ’)-n e 


(2m e ) 


Vi 


(42) 


3/2 „ - 3f° _ 3f° 

e 3/2 V x u ^ - W 


COLL 


*In McKenzie’s derivation the expression for A g k,k 1 should be 
replaced by that given by Eq. (13). 
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and 



(n e f ') + V x (n e f°) - n e eE j- -n e m e ^ 7 - 


3f° 


Dt 


Be 


Du 

Dt 


3f° 

Be 


BT 

Bt 


'coll 


(43) 


where f is normalized such that 

r 00 ^ 

J e h f°(e,x ,t)de = 1 . 


(44) 


With this normalization, if the distribution were Maxwellian, at temper- 
ature T (*,t), under the present approximation the distribution function 

would bl given by 


fmax ( e i x >t) 


2 ]_ 

(«T e ) 3/2 


exp(-e/»cTe) . 


(45) 


Note that from Eqs. (42) and (43) the conservation equations of electron 
energy and momentum can be derived. To obtain the electron energy equation 
multiply Eq. (42) by (±_) 1/ ' 2 ede and integrate over all electron energy. 

With the electron energy equation written in such a form it is then 
possible to obtain the fraction of power transferred by the electrons 
to the various modes of the gas, in particular, to the vibrational levels. 
In this way the parameters in the problem can be adjusted to maximize 
the amount of power transfered to the vibrational levels and concurrently 
the laser efficiency. Such an analysis has been carried out for a 
static discharge by Nighan 27 . To obtain the momentum equation multiply 
Eq. (43) by 2/3 e 3 ' <le and integrate over all electron energy. 


Collision Terms 

The remaining expressions required are those for the collision terms. 
Considerable simplifications arise in these terms because the electron 
mass is much smaller than that of the heavy particles, causing the heavy 
particle speed to be much smaller than that of the electrons even for 
high supersonic speeds. Thus, the collision operators can be approximated 
by the operators occuring in a plasma of motionless heavy particles an 
moving electrons. For the energy transfer mechanisms included in this 
model the collision operator consists of a sum of four contributions 
accounting for elastic-heavy particle collisions, 8f / 9t | e h, ionization 


and three-body recombination collisions, ^i/gtlioh. 


charged particle 


collisions 9f/ 3t i_ e and most importantly, inelastic and superelastic 

excitation and de-excitation collisions / 3t | ex ' Expressions for the 

collision terms for these types of collisions emplo^n^the small electron 

mass approximation have been developed in the past an 

those appropriate to the present case will be given. For the f operators 

we have: 
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(46) 


at 


eh 


= "* ? ST A { e2 ">Q“ h < E > [f°W + kT } , 


3f° 

at 


s f r e-ej 

n e Z1 X) I e C Sk nj n e J Ss X (e,e b )f°(e b ) • f°(e-e b -e^)de b 

s k=0 l O 

r 00 

-en k f°(e)S kx (e) + 2.n k J ejf 0 (e i )S s kx (e i> e)de i 

o 

00 

-2n‘n e C Sk f°(e) f e i f°(e i -e-6 kx )S kx (e i ,e)de i 


(47) 


3f° 

at 


n e 2^j3n e (»cTe) 2 Qr 


af°(e) n o + i° n 

be (h+J-i) 


1 [1; 


f° (€) + je 


3f° 

at 


N s— 1 N s [ 

= n e Z 2 2 | (e+e kc )S kc (e+e k8 )n k 


s k-o fi=k+l 


rf°(e+£, k *) - 44 **wl + " 

L g s 4 

-HI • 


k eS k8 (e) 


(48) 


(49) 


r gs ( h 2 \ 3/2 

% 2g x \ 27rrrie/ 


vX 
2 ’ 


where f°(e,x, t) is written as f°(e) for convenience. 

Here the sum in Eq. (46) includes all heavy particle species, and 
e =£ kX +e +£ t ^ e energy of the incident ionizing electron. The 
i s b 

following definitions also apply: 

= energy-dependent electron-heavy-particle momentum 
transfer cross section, 

= energy average electron-ion cross section, 

= 6n(e 2 /12lIe o kTe) 2 £nA, 


qf(e) 


,ei 
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degeneracy of level k of species s ; unity for vibrational 
levels, 


S g (e,e^) = Cross section for ionization of a particle in the kth 
level by an electron of energy e which results in an 
electron of energy e, . It is related to the conven- 

b kx 

tional ionization cross section S (e) as follows: 

e-Cj x S 1 

S s k "(e) = / S S kk (e,e b )de b , 

o 

lp--L 

e p 2 6 

Jp=-V de • 

e p/2 { 

Note that these expressions are written for a general level k. In other 
words, with the proper interpretation of the cross sections, the equations 
apply for either rotational, vibrational or electronic level excitations. 
However, to apply them to situations which involve all three processes 
simultaneously, it must be assumed that the cross section for each process 
is independent of the other. In this way each process can be described 
by equations written in the above form and each expression can then be 
added to obtain the total contribution. 

To determine the expressions for the? collision operator. Refs. 44 
and 46 show that only the elastic collision term is important. This is a 
consequence of the fact that the electron^heavy-particle momentum 
transfer cross section is larger than any of the other cross sections. 

With this we find that the f' operator can be written in the following 
simple form: 

= ~ n e [fn s Q® h (e)] f 1 (e,x ,t) , ( 50 ) 

COLL 


where the sum includes all heavy particles. 

The necessary expressions for the collision terms have now been 
developed for the Boltzmann equation and, by making certain order of magnitude 
arguments applicable to the present model, certain simplifications are . 
possible. By a comparison of the magnitudes of the various cross sections 
(the sources of these cross sections will be given latter), we find, for 
the weakly ionized gas assumed here, with m e <<m g , and for the average 

electron energies involved, that: a) the elastic collision term Eq. (46) 

can be neglected, b) the first: and fourth terms in Eq. (47), which 
represent recombination terms can be neglected, and c) the charged- 
particle collision term Eq. (48) can be neglected. 
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Simplified Equations 


We seek to further simplify Eqs . (42) and (43) by examining the 
various terms on the lef-hand side of the equations, where we now 
specialize to a steady one— dimensional case. To facilitate this we 
introduce the following non-dimensional variables: 


— u — 

= _P_ j 

e — 

u ■ 

Po’ 

2 5 ^ 
VlTTiQ Cq 

vks _ X* 2 (e) 
Xs ' X R 

, x eh = 

X eh (e) Yk x _ 

Xr '* Xs 


X^(ei,e) 
X R (VimeCo) 


F = 


F 

Po (Vim e Cq ) -3/2 


where 


X k8 (e) 


1 

S ke (e)p/m 


is an inelastic collision 


mean free path. 


A eh (e) 


?n s Qf 


is an elastic collision 


mean free path, and F = px f. The above normalization is chosen such 
that the barred variables are of order unity with the exception of F, 
where F’«F°. With this normalization Eqs. (42) and (43) become, after 
substitution of the simplified collision terms, the following: 



(51) 


and 

/ Xr \ / Uo\_u_ d F 1 + / Xr\ d F 0 
V L / We* dy \ L / dy 

= - F 1 (e)/ A eh , 

where 



- 2 



-du 3F° 
u — — — — 

d y de 


(52) 


dF° 

dt 


. (e+e s 8 ) k 
COLL = ? “ 5 « X ‘ 


r> o / i 

F (e+e s ) 


F° (e)l + _ X ^ 6 
e «x k J X k8 (i) 


g s 

, ex k F°(i) 

+ 2 2 \ — + x, 

S 


X^ x (e) 


F°(e-e ks )^ - F°(e) 

gs Xs 


o X kk (e i; e) 


(53) 
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Three parameters appear in these equations, namely 


2eEA R 

LH e Cq 


and 




Upon estimating the magnitude of these quantities we find that for the 
average electron energies involved in this problem (1 ev) u q ^ <<1 > 

o 

and the ratio 2 e EA R 


is of order unity. The ratio A R /L, which is the 


m c 

Knudsen number, fs°in general less than one and we will exploit this by 
expanding the distribution function F in a series in A R /L=6. This 

expansion is given by: 


F° -> F +SF + . . 
o 1 

where 6<<1. 


(54) 


Thus, taking into account the above order of magnitudes, substituting 
Eq. (52) into Eq. (51) to obtain an equation for F alone, using the 
expansion (54) in the resulting equation, retaining only the first two 
terms of the expansion, we obtain the following two coupled equations for 
F and F (dropping bars for convenience) 

^ i \ ^ C I 

(55) 


-T E « 


_9_ 

1 c > eh 

aF 0 1 

9F 0 

de 

r 

96 | 

9t 

a 

L^eh 

9F, \ 

9F i 

9e 

r 

de 

" 9t 


'coll 


2 ./Uo\ , 3/2 du 9 fo _ / Uo \ M 
3 \C 0 J dx 3e \ c o/ 


'coll 
9F 0 


_ J_p2 9 

9x 3 R £ 9x 


l eh 


9F 0 

9e 


(56) 


J_ F 2 2L 

3 Er 9e 


where 


3Fo 

at 


eA eh 

and 
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9F 0 

dx 

3Fi 

~aT 


'coll 


»v ^ 1 r ' 

for F 1 , to consistent order in 6 and u o/c ^ is given by: 


with $he appropriate substitution, and 

in 

9F 0 


_ 2eEA R 
Er = ~ 


F 1 = - A 


-Er 


3e 


are given by Eq. (53) 

Finally, the equation 
(57) 


m e Co 


+ 5 


dF 0 3Fi 

3x R de 


Equations (54)- (57) give the solution to the free-electron Boltzmann 
equation provided the Knudsen number is small. This confines the 
region of applicability of these equations to situations where gradients 
in the flow are not large. 

We see from these equations that the effects of gradients and flow 
velocity are demoted to second order. In fact, as one might expect for 
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such a situation the equation for F , which gives the largest contribution 
to the distribution function, is the same as that for a stationary gas 
(see Nighan^ ). The difference here is that Eq. (55) is to be interpreted 
as being applied locally at each point in the discharge. In a gross sense, 
one can then say that to first order the presence of the discharge alters 
the flow, but the flow does not alter the electron distribution in the 
discharge. 

To solve these equations the cross sections for electron-vibrational, 
electronic, and ionization collisions important to this problem must be 
known. A search of the literature gives many of these cross sections 
for various gases that might be employed and certain of these references 
are listed below: 

CO: vibrational excitation (ref. 55), electronic excitation (ref. 56, 

61), momentum transfer (ref. 56), ionization (ref. 62). 

N : vibrational excitation (ref. 55, 57, 59) electronic excitation 

and momentum transfer (ref. 57), ionization (ref. 62). 

He, Ar, Xe: electronic, ionization and momentum transfer (ref. 58, 

60, 62). 

For the crucial electron-vibrational collisions in CO only cross 
sections for excitation from ground state molecules are known (ref. 55). 
Thus, when including electron collisions with excited state mg^ecules 
the cross sections must be estimated (see Abraham and Fisher for one 
scheme of estimation) . 

Finally, in order to solve Eq. (55) (and(58)) which is a second-order 
ordinary differential-difference equation, a numerical technique is 
required. Several approaches have been used in the past 4 ’ ’ each 

with its own virtues and the interested reader can refer to the references 
cited for more information. 

With a solution procedure for the free-electron Boltzmann equation 
chosen, the remaining question is the approach to be used to obtain the 
numerical solution of the coupled set of rate, gasdynamic, and electron^ 
kinetic equations. A possible scheme is to employ the method of Bailey 
who obtained solutions for nozzle flows with chemical nonequilibrium. This 
scheme linearizes the gasdynamic and rate equations and then employs an 
implicit integration scheme to obtain the solution of these linearized 
equations. This method has been used successfully by McKenzie in his 
gasdynamic laser studies. 

When using such a scheme the electron properties need not be considered 
as independent variables such as T, u, etc., and thus they become 
auxiliary quantities to be determined at each point. This implies that 
the free-electron Boltzmann equation solution also becomes an auxiliary 
calculation,, and this can greatly reduce the complexity of the numerical 
calculations . 
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